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We consider an integrable generalization of the nonlinear Schrodinger (NLS) 
equation that was recently derived by one of the authors using bi-Hamiltonian 
methods. This equation is related to the NLS equation in the same way that the 
Camassa Holm equation is related to the KdV equation. In this paper we: (a) 
Use the bi-Hamiltonian structure to write down the first few conservation laws, 
(b) Derive a Lax pair, (c) Use the Lax pair to solve the initial value problem, 
(d) Analyze solitons. 

AMS Subject Classification (2000): 35Q55, 37K15. 

Keywords: Integrable system, inverse spectral theory, Riemann-Hilbcrt problem, solitons. 

1 Introduction 

The celebrated nonlinear Schodinger (NLS) equation appears in a wide range of phys- 
ical circumstances from water waves to nonlinear optics. This is a consequence of the 
fact that it can be derived using physically meaningful asymptotic techniques from 
large classes of nonlinear PDE's [2 . Indeed, the NLS equation describes the ampli- 
tude A(£, t) (£ = ex, t = e 2 t) of the modulation of a monochromatic wave of a large 
class of dispersive nonlinear equations. 

NLS belongs to a distinctive class of nonlinear equations called integrable. Inte- 
grable equations possess several special features including the existence of a Lax pair, 
as well as a bi-Hamiltonian formulation. Actually, the latter feature provides a simple 
algorithmic approach of constructing integrable equations. Indeed, the following re- 
sult was derived independently in [5] and [T21 E3] : Suppose that the operator 9\ + A#2 
is Hamiltonian for all values of the parameter A (i.e. this operator is skew-symmetric 
and it satisfies an appropriate Jacobi identity). Then, each member of the hierarchy 
of equations 
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n G Z, Cj constant 



possesses infinitely many conserved quantities. Equation (1.1) with c\ = —1, N = 1, 
and 

(1.2) 6q = <9, 6 2 = d + f3d 3 + -a(qd + dq), a, f3 constants, 

3 

yields the KdV equation 

qt + q x + /3<2W + aqq x = 0. 



Similarly, equation (1.1) with c\ = —1, N = 1, and 

6q = <9 + z^(9 3 , #2 as in (1.2), ^constant, 



yields, with q = u + vu xx , the Camassa-Holm equation 

(1.3) u t + Ma + vu txx + /3n a . a . x + auu x + ^au(uu xxx + 2u x u xx ) = 0. 

The KdV equation provides an approximation to two-dimensional inviscid, irro- 
tational water waves above a horizontal bottom y = —ho under the assumptions that 
the parameters a = aj ho and f3 = h^/l 2 are sufficiently small, where a and I are 
typical values of the amplitude and the wavelength of the waves. In these circum- 
stances, one obtains the KdV if 0(a) = 0(/3) and one neglects terms of 0(a 2 ). It 
is shown in |1U] that if one retains terms of 0(a 2 ) (but neglects terms of 0(a 3 )), 



one obtains the Camassa-Holm equation. Although equation (1.3) was introduced by 
Fuchssteiner and one of the authors in 1981 |llj . it began to attract a lot of interest 
(see for example |U SI El ITU H71ll9j) only after its rederivation by Camassa and Holm 

It was noted earlier that it is possible to derive mathematically the CH equation by 
utilizing the two Hamiltonian operators associated with the KdV equation. Similarly, 
by utilizing the two Hamiltonian operators associated with the NLS equation, it is 
possible to derive |6] the following generalization of the NLS equation 

(1.4) iu t — vutx + ^yu xx + a\u\ 2 (u + ivu x ) = 0, a = ±1, 

depending on two real parameters 7 and v. In the same way that the CH equation 



provides a better approximation of water waves, equation (1.4) appears in the same 



circumstances as NLS provided that one retains terms of the next asymptotic order 



|6j. In this paper we: (a) Review the derivation of (1.4) and write down the first 
few conservation laws using the bi-Hamiltonian structure, (b) Derive a Lax pair of 
equation ( |1.4[ ). (c) Use the Lax pair to solve the initial value problem, (d) Analyze 
solitons. 

2 Derivation and bi-Hamiltonian structure 

Consider the NLS equation 

(2.1) iqt + jqxx + crq\q\ 2 =0, a = ±1, 
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where 7 is a real parameter. This equation is a particular reduction of the following 
bi-Hamiltonian system: Define the compatible pair of Hamiltonian operators 



the system 



1 \ / jd x + qd x V -qd x 1 q 

cr 3 = [ and 92= ' „ , x a l 1 

\0 -1/ V -rd x r ^yd x + rd x l q 



reduces to equation (2.1 ) when r = aq. If we replace 03 by 0\ = 03 + ivld x where / 



is the 2x2 identity matrix, then the equation 



together with the definition q = u + iuu x , gives (1.4). 
Now let 

(0 1 

ai = (,i 

and define two Poisson brackets for the functionals F[q, r] and G[q, r] by 

{F,G} V = /(gradF) T PgradGdx, 



{F, G} £ = J (grad F) T £grad Gdx, 



where the superscript T denotes the transpose and the operators T> and £ are defined 
by 

-qd x l q >yd x + qd x V 



2? = 2 ai 



jd x + rc^ : g — rd x 1 r 



and 



-lV\(Tl 





i + ud x 



-i + vd x 




The gradient of a functional F[q, r] is defined by 



g radF= , 

provided that there exist functions ^ and ^ such that 
d 



de 



F[q + eipx,r + e^ 2 ] 



<E=0 



VV 1 + ^ 2) ' 



for any smooth functions <*p\ and y>2- Integration by parts shows that the operators 
T> and £ are skew-symmetric with respect to the bracket 



(2.2) 



(hgi + hd2)dx. 
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If we introduce u and v by the relations 
(2.3) q = u + iisu x , 



r = v — iWr, 



we find that equation (1.4) has the bi-Hamiltonian formulation 

= Pgrad Hi = £grad H 2 , 



where the Hamiltonians Hi and H2 are given by 



H\ = i I u x rdx = —i v x qdx 



and 



Ho 



/( 



^vuxx + -quv z )dx 



/( 



1 



'yuvxx + -ru v )dx. 



Indeed, if we write equation (|1.4|) as the system 

lUxx + quv 
-IVxx ~ ruv 



then using 



grad Hi 



-%v x 
iu x 



we immediately find the formulation in terms of H\. Moreover, from (2.3) and the 



definition (2.2) of the bracket, we infer that 



grad H 2 



Hence, using 



SH 2 
5u 

6H 2 
8v 



l-iud x )- 1 SH 2 /Su 
1 + ivd x )- l 5H 2 /5v 



-jv xx - ruv 
~iu xx - quv 



we find the formulation in terms of H 2 . A hierarchy of conserved quantities H n can 
be obtained by means of the recursive relations 



We find that 



is satisfied for 



VgradH n = Sgvad H n+ i, n £ Z. 



Dgrad Hq = £grad Hi 



Hq = — I rqdx. 
7 J 

The first few conservation laws and associated Hamiltonian equations when q = or 
are summarized in Table [TJ 
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Table 1: Conservation laws of equation (1.4) 



n 


H n 


q t = £ grad H n 


2 


- J (yvu xx + \quv 2 ) dx 


qt = iju xx + iaq\u\ 2 


1 


i J u x r 


qt = q x 





i J rqdx 


qt = -i(l+iud x )q 



3 A Lax pair 

Replacing u(x, t) by u(—x, t) if necessary, we can assume that v has the same sign as 
7 in equation (|1.4[). Let a = 7/1/ > and (5 = 1/za Then the transformation 



-<r, 



converts (1.4) into 
(3.1) 



utx + ol (3 u — 2ia(3u x — au xx + oiafi \u\ u x = 0, 



a = ±1. 



Thus, without loss of generality, in the remainder of this paper we will concentrate 



on equation (3.1 ). This equation is the condition of compatibility of the system 
ip x + i( 2 a 3 ifj = (U x i/j, 



(3.2) 



aCUx + ic^^(l U _ U 2 



where ip(x, t) is a 2 x 2-matrix valued function and 



(3.3) 



U 



u 
i? 



V = V<x[( 



2C 



u = au. 



Notice that the x-part of (3.2) has the same form as the rc-part of the Lax pair for 



the derivative NLS equation derived in the pioneering work of Kaup and Newell [16 . 



The fact that the x-part of (3.2 ) contains only the x-derivative of U is related to the 



absence of a i^-term in the gauge transformed equation (3.1). 



4 The solution of the Cauchy problem 

In this section we will use spectral analysis to express the solution of the Cauchy 



problem of (3.1 ) in terms of the solution of an appropriate Riemann-Hilbert problem. 
Letting 

(4.1) il>(x,t,Q = e^-« v -^' t ^ ( ^' t ^^M(^,t s C)^(s^)e^ a ^* )ws , 

where 03 acts on a 2 x 2 matrix A by a%A = a^Aa^ 1 and the diagonal matrix D is 
defined by 



(4.2) 



D(x,t) 



e 2 



y™v x (x>,t)u x (x>,t)dx'a 3 
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the Lax pair (3.2 1 becomes 



(4.3) 

where 

(4.4) V l 

(4.5) V 2 



fit +ir/ 2 [(T3,/x] = V 2 n, 



e -U-o v ^ dx ' a 3 



e -U- a o V * U * dx ' a 3 



-v x u x a 3 



1. 



a<U x + -^-a 3 l-U - I 1 



■ a I n2 \ 

i-^{u x v x - (3 uv)a 3 



In order to derive (4.3) from (3.2) we use the conservation law 
(4.6) {u x v x ) t = a{u x v x - (3 2 uv) x 

to carry out the differentiation of D and of the exponential prefactor on the right-hand 



side of (4.1 ) with respect to t. 



Before using the spectral analysis of (4.3) to derive a nonlinear Fourier transform 



pair, we note that the transformation (4.1) can be motivated as follows: If U x tends 
to zero at infinity, then there exists a solution of (3.2) which tends to exp(— i(£ 2 x + 



7j < t)a 3 ). This suggests the transformation ip = \I/exp(— i(£ x + rj t)a 3 ) and then ^ 
solves 

(4.7) 



* x + iC 2 [a 3 ,V] = (U X V, 



aCU x + ^fa 3 UU-U 



We seek a solution of this equation in the form 
where D, $1, ^2 are independent of (. Substituting the above expansion in the x-part 



00, 



of (4.7) it follows from the 0(C) terms that D is a diagonal matrix. Furthermore, 



one finds the following equations for the 0(C) and the diagonal part of the O(l) terms 



-U x Da 3 , 



(4.8) 0(():i[a 3 ,y 1 ]=U x D, i.e. 

where denotes the off-diagonal part of ^1; 

0{l):D x = U x ^ { °\ i.e. D x = % -u x v x a 3 D. 



This leads to the definition (4.2 ) of D. In fact, a similar analysis of the i-part of (4.7 ) 
reveals that Dt = ^t(u x v x — 2 uv)a 3 D, which, in view of the conservation law (4.6), 



is compatible with (4.2). Therefore, up to a multiplicative constant, the asymptotic 



behavior of Vl/ is given by D as defined in (4.2). In order to formulate a Riemann- 



Hilbert problem for the solution of the inverse spectral problem, we seek solutions of 
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the spectral problem which approach the 2x2 identity matrix I as £ — * oo. This 



suggests introducing /j, by ip = /iDe 



-i(<; 2 x+ v 2 t)<73 



3 , so that 



(4.9) 



n = i + o 



oo. 



The additional exponential factor in (4.1), which does not affect (4.9), is convenient 



because it implies that two different solutions of (4.3) are related by multiplication 



by a matrix which is independent of U x (see equation (4.14) below). 



4.1 A nonlinear Fourier transform 



We now consider the spectral analysis of the x-part of (4.3). Since the analysis will 



take place at a fixed time the i-dependence will be suppressed. Define two solutions 



and [i2 of the x-part of (|4.3|) by 
(4.10) 

(4.11) H2(x,0 = I 



Ml 0,0 = 1 + I 

oo 
oo 



iC, 2 (x'~x)a 3 



J( 2 (x'-x)& 3 



(Vini)(x',()dx', 
(V lf i 2 )(x',()dx'. 



The second columns of these matrix equations involve exp[2i( 2 (x' — x)]. It follows 
that the second column vectors of \x\ and n% are bounded and analytic for £ in D- 
and .D+, respectively, where 

D + = {CI arg C G (0, 7r/2) U (tt, 3tt/2)} 

is the union of the first and third quadrants, while 

£>_ = {C| arg C G (vr/2, vr) U (3vr/2, 2^)} 

is the union of the second and fourth quadrants. We will denote these vectors with 
superscripts (— ) and (+) to indicate these boundedness properties. Similar conditions 
are valid for the first column vectors, hence 



Mi = } ) , M2 

Note that fj,(x,Q = fj,j(x,Q, j = 1,2, satisfies the symmetry relations 
(4.12) 
as well as 



a4 ^^4 



MuO, = M22O, cO, fjvi(x, C) = M12O, O"0, 



(4.13) 



M11O) -C) = Mnfo 0, M12O, -C) = -M12O, C)) 

M21O, -C) = C)> M22O, -C) = M220, 0- 
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4.1.1 The function s(() 

Any two solutions fi and jl of ( |4.3[ ) are related by an equation of the form 
(4.14) Mx,C) = /2(x,C)e-^ 3 iVo(C), 

where iVo(C) is a 2 x 2 matrix independent of x. Indeed, let ip and ip be the solutions 
of 



(4.15) 

related to fi and jl as 



li{xX)De- K2xa \ 



Then, since the first and second columns of a solution to (4.15) satisfy the same 
equation, there exists a 2 x 2 matrix iVi(0 independent of x such that 

(4.16) ^(x,C) = feC)JVi(C). 

Comparing equations (4.14) and ( |4.16 ) we find 

N (0 

Thus we can define the spectral function s(0 by 



e - 1 JTL WxWa^o-a iVi (C) . 



(4.17) 

Evaluation at x 
(4.18) 



/i 2 (x,C) = w(^C)e- jC x&3 s(0, 
— oo gives 



ImC 2 



0. 



s(C) = I 



^ 2xa3 (V lf i2)(x,()dx, ImC 2 



0. 



Moreover, it follows from (4.3) that the determinant of /i is independent of x, so that 
evaluation of det /i,-, j = 1,2, at x = ±oo shows that 

(4.19) det^ = l, J = 1,2. 

In particular, 

dets(C) = I- 



From (4.12) and (4.13), we deduce that 



Sll(C) = 322(0, s 2l(C) = 0-512(0- 

This analysis suggests the following notation for s: 



«(C) 



^6(0 «(0, 



From the explicit expression (4.18) for s(0 and the fact that the second column of 
[12 is defined and analytic in we find that a(0 has an analytic continuation to 
all of D + . The symmetry (4.13) implies that a(0 is an even function of (, whereas 
6(0 is an odd function of £• 
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4.1.2 Residue conditions 



Since a(£) is an even function, each zero (j of a(£) is accompanied by another zero 
at — Q. We assume that a(£) has 2N simple zeros {Cj}|=! C D + such that {Cj}jLi 
belong to the first quadrant and Cj+N = — £j) j — ,N. 
The second column of equation (4.17) is 



(4.20) 



/4 +) 



a/4 ^ + b^'e 



( £ 1U iR. 



Applying det I fJ^\ -J to (4.20) and recalling (4.19), we find 

det (/z£ +) (s,C),^ +) (z,C)) = «(0, C G £>+, 

where we have used that both sides are well-defined and analytic in D + to extend 
the above relation to all of D + . Hence, if a(£j) = 0, then /J^\x, Q) and (sc, £,-) 
are linearly dependent vectors for each x. It follows that the 2x2 matrix 



is a solution of the x-part of (3.2) with linearly dependent column vectors. We 



conclude that there exists a constant bj such that 

= b j ^ + \x,( j ), xeR. 

This implies 

(4.21) fi[ + \x, (j) = bje 1 ^ v * u * dx e 2i $ x ij£\x, Q, x G R. 

Recalling the symmetries ( |4.12 )-(4.13), the complex conjugate of (4.21) is 



(— )/ j- \ 7 — j f°° v x u x dx — 1i(?x (— )/ > \ 



x e JR. 



Consequently, the residues of /i^ (x, C)/a(£) and • ) (x,C)/o(C) a t Cj anci Cj are 

,^i +) (x,o _/4 + Wi) 



Res- 
Cj o(C) 

.(-)/ 



HQ) 



<* a(C) 



where d(C) = and 



1 i f°° v x u x dx 



a(Cj) 



In fact, there are only N independent constants Cj. Indeed, using the symmetries 
(4.13) it follows from (4.21 ) that frj+jv = — bj, j = 1, . . . ,N. The fact that d is an odd 
function implies the relations Cj+N = Cj, j = 1, 
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4.1.3 The inverse problem 



The inverse problem invo 



ves reconstructing U from the spectral functions fj,j(x,Q, 



j = 1,2. We obtained in (|4.8|) that = ~U x Da 3 whenever 



^1 #9 / 1 \ 



oo, 



is a solution of (4.7). This implies that 

(4.22) u^x) = 2im{x)e i S-°° VxUxdx \ 

where 

ttjC 1 ) m ( 2 ) / i 



C ^ oo, 



is the corresponding solution of (4.3) and we write m{x) for m^(x). Prom equation 
(4.22) and its complex conjugate, wg obtain Vux^x 

= 4a\m\ 2 . Thus, u x is expressed in 



terms of m by 
(4.23) 



u x (x) = 2im(x)e 



Aai [ x \m\ 2 dx' 



4.1.4 The Riemann-Hilbert problem 



Equation (4.17) can be rewritten in a form expressing the jump condition of a 2 x 2 



RH problem. Algebraic manipulations give 

M_0,C) = M+(x,QJ(x,0, CeKU«, 
where the matrices M_, M+, J are defined by 



(4.24) M 4 



fi 2 



(-) Mi 



a(C)/ 



(4.25) 



J(x,C) 



fc(C) e -2»c 2 ^' 

a(C) 



C 6 M U iM. 



a(C)a(C) 



The contour for this RH problem is depicted in Figure [TJ We summarize our discussion 
in the following theorem, in which we use the notation [A]\ ([A] 2) for the first (second) 
column of a 2 x 2 matrix A. 

Theorem 4.1 1. Given u x (x) with sufficient smoothness and decay, define the 
map 

§:{u x (x)}^{a(0,b(0,Cj} 
by the following equations: 



(4.26) 



6(0 
a(C) 



»(C)]2, s(C)= lim ''"''' !; V'2(.r.C). CGKUfl 
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J 



D_ 






D_ 



Figure 1 Illustration of the contour and the four sections of the Riemann-Hilbert 



problem (4-24) in the complex (-plane. 



and 

(4.27) Mx,Ci)] 2 a(0)C,- = Mx,0)]ie- 2 ^, j = l,...,2N. 

The functions fj,i(x,Q and fi2(x,Q are the unique solutions of the Volterra 
linear integral equations (^7o)-(^T7p ; and we assume that a(C) has 2N simple 
zeros {Cj}j=! C D + ordered so that {CjljLi belong to the first quadrant and 
Q+N = ~Cj, j=l,...,N. 

The spectral data {a(C), 6(C), Cj} have the following properties: 

• a(C) is defined for ( G D + and analytic in D + . 

• 6(C) is defined for Im C 2 = 0. 

• a(CMC)-^(C)RC) = l, ImC 2 = 0. 

• o(C) = l+(j), C^oo, (£D+- 

• 6(C) = (J) , C^oo, ImC 2 = 0. 

• a(C) 6(C) satisfy the symmetry relations 

(4.28) o(-C) = o(C), CG^+, 6(-C) = -6(C), ImC 2 = 0. 

• C j+N = Cj, j = 1,...,N. 

2. The map Q : {o(C), 6(C), Cj} i— > {u x (x)} ; inverse to S, is defined by 

(4.29) u^x) = 2im(x)e 4CT ^-- |m|2d:r ', m(s) = lhn (CM(x,C))i2, 

>oo 

where M(x, C) «s i/ie unique solution of the following RH problem: 
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M(x,C) 



M + (x,() (ED + , 
M_(x,C) 

is a sectionally meromorphic function. 

M-(x, C) = M + (x, ()J(x, C) for ( £ R U iR, where J is defined in H4.2ty . 
M{x, Q has the asymptotic behavior 



(4.30) 



M(x,C) =1 + 



C 



oo. 



The first column of M + has simple poles at Q = j = 1, 
the second column of M_ has simple poles at £ = j = 1, . 
associated residues are given by 



,2N, and 
,2N. The 



(4.31) 
(4.32) 



Res[M(x, C)]i =C,e 2 ^[M(x, C;)] 2 , 



Res[M{x, C)] 2 =C J e- 2 ^ :c [M(x, (.•)] 



j = l,...,27V, 
j = l,...,2iV. 



5. We Ziaue 



3-1 



In the case when a(£) has no zeros, the unique solvability of the Riemann-Hilbert 
problem in Theorem |4.1| is a consequence of the following vanishing lemma, while if 
a(Q has zeros the singular RH problem can be mapped to a regular one following the 
approach of [9] . 



Lemma 4.2 (Vanishing lemma) The Riemann-Hilbert problem in Theorem 
with the vanishing boundary condition 



M(x,()^0 as C -> °o, 



has only the zero solution. 



Proof. Assume that M(x,Q is a solution of the RH problem in Theorem 4.1 such 
that M±(x,Q —>■ as £ — > oo. Let A* denote the complex conjugate transpose of a 
matrix A. Define 

H+(() = M + (C)Mt(-a(), (£D + , 
F_(C) = M_(C)M|(-aC), CeL>-, 

where we have suppressed the x dependence. The functions -H+(C) an d -ff-(C) are 
analytic in the interior of D + and D-, respectively. By the symmetry relations (4.28) 
we infer that 

jt(-aC) = J(C), Ce^uiR. 

Since 



H + (() = M + (()j\-o-()Ml(-o-(), H-(C) = Af+(C) J(C)M|(-aC), C G K U ill 
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it follows that H+(C) = H -(0 on tU«. Therefore, H + (() and #-(0 define an 
entire function vanishing at infinity, and so -ff+(0 an d H-(C) are identically zero. 

Let us first consider the case a = 1. In this case J(is) is a hermitian matrix 
with unit determinant and (11) entry 1 for any s. Hence, J(is) is a positive definite 
matrix. Since H_(is) vanishes identically for s £ 1, we find 

M + (is)J(is)M ] + (is) =0, s£l, 

and so M + (is) = for s£i By analytic continuation, M + and M_ vanish identi- 
cally. This proves the lemma when a = 1. 

If a = — 1 , we note that J(s) is a hermitian positive definite matrix for any s£t. 
Since -ff-(s) vanishes identically for s £ 1, we find 

M+(s)J(s)Ml(s) = 0, seR, 

and so M + (s) = for s £ M. Again it follows that M + and M_ vanish identically, 
which proves the lemma when a = — 1. □ 



Remark 4.3 It follows from the symmetries (4.12 )-(4.13) that the solution M(x,0 
of the Riemann-Hilbert problem in Theorem |4.1| respects the symmetries 



(4.33) 



Mn(x,C) = M 22 (x,a(), M 2 i{x,() = M 12 {x,a(), 
M u (x, -0 = M n (x, 0, M 12 (x, -0 = -M 12 (x, 0. 
M 21 (x, -0 = -M 21 (x, 0, M 22 (x, -0 = M 22 (x, 0. 



Moreover, if one enforces (4.33), then only one of the two residue conditions (4.31)- 
(4.32) needs to be verified since the other condition is a consequence of symmetry. 



4.2 Time evolution 



Suppose u(x,t) is a solution of equation (3.1). At each time t, we define the spectral 
data a((,t), b((,t), and {C^t)}^ according to \A.2%\ and \A.27\ , i.e. 



(4.34) 



KC,t) 

o(C,t) 



[s(C,t)} 2 , s(C,t)= lim e^ xa ^ 2 (x,t,C), (GRUtl 



and 
(4.35) 



[A*2(x,t J 0(t))] 2 A(0(*).OCi(t) = MM,0(t))]ie- 2 ^ (t)x , j = l,...,2iV. 



Proposition 4.4 The evolution of the spectral data is given by 
(4.36) 

o(C,*) = o(C,0) ) 6(C,t)=6(C,0)e- 2 ^, 0(*) = Ci(0), ^-(t) = C^Oje 2 ^*, 



where rj was defined in (3.3) and rjj = \fa (0 — gf" 
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Proof. From the i-part of (4.3), we infer that 

(4.37) (e i ^ x+ ^ a ^ 2 {x : tX)) t = e i(cWi) " 3 (^ 2 )(x,t,C). 

Assuming that u has sufficient decay at infinity, we have V2 — > as x — ► ±00. 
Evaluation of (4.37) at x — > —00 gives 



0, 



C £ RU iR. 



The second column of this equation gives the time evolution of a and b as stated in 
(4.36). Since at = 0, the eigenvalues O are time- independent. Moreover, differentiat- 



ing the first row of (4.35) with respect to t and evaluating the resulting equation at 
x — ► —00, we find 

,2/ 



C jt = 2iTfC j , j = l,...,2N. 



This completes the proof. 



□ 



4.3 The Cauchy problem 

The Cauchy problem for equation ( |1.4[ ) with initial data uq(x) with sufficient smooth- 
ness and decay can now be solved using only linear operations as follows: 

1. Define fi\(x,Q and /^(z, C) m terms of uq x as the unique solutions of the 



Volterra linear integral equations (4.10 )-(4.11 ), and use them to obtain the 
spectral data {a(C, 0), b((, 0), Cj(0)} at time t = according to (|4.26|-(|4.27[). 



2. Obtain the spectral data {a(£,t),b(£,t),Cj(t)} at a later time t from equation 

(pi. 



3. Solve the RH problem in Theorem 4.1 with the jump matrix J(x, t, C) defined in 



terms of a((, t) and 6(£, t) by (4.25 ) and with the residue conditions (4.31 )-(4.32 ) 
defined in terms of Cj(t). 



4. Use the solution M(x, t, Q to obtain u x (x,t) according to (4.29). 

5. Obtain u(x,t) by integration. 



Observe that (3.1) is an evolution equation in u x and that any solution u(x,t) 

u(x,t) + h{t) for an arbitrary function h(t). The 
00 removes this non-uniqueness. 



is undetermined up to u(x, t) 
requirement that u goes to zero as \x 



5 Solitons 

The solitons correspond to spectral data {a(C)> HC)j Cj} for which 6(C) vanishes iden- 



tically. In this case the jump matrix J in (4.25) is the identity matrix and the 



Riemann-Hilbert problem of Theorem 4.1 consists of finding a meromorphic function 
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M(x, £) satisfying (4.30) as well as the residue conditions (4.31 )-(4.32). From (4.30) 

2iV „ 2iC 2 x\- hlT( > m 

Cj-e ^ [M(x,(j)\2 



and (4.31 ) we get 
(5.1) 



[M(x,0}i=(l)+J2 



3=1 



c-o 



If we impose the symmetries (4.33), equation (|5.1|) can be written as 

2N 

Evaluation at a^k yields 



v r^s- /M 12 (x,Q) 



(5.2) 



M 22 (x,q) 
M 12 (x,Cfc) 



2A ! 



v (> a 'V /M 12 (x,C, 



fe = 1,...,2N. 



Solving this algebraic system for Mi 2 (x, Q), M 22 (x,(j)i j = 1,...,2N, and substi- 
tuting the solution into (5.1 ) yields an explicit expression for [M(x, C)]i- If we define 
[M(x,C)]2 by (4.33), the second residue condition (4.32) is a consequence of symme- 
try. This solves the Riemann-Hilbert problem. The potential u x {x) can be found by 
noticing that (5.1) implies 

2N 



M 12 (x, ()=<rJ2 C 1 e- 2i ^ x M 22 (x, Q)~ + O 

3=1 ^ 



Therefore, by (4.29), 



(5.3) u x (x) = 2im(x)e 



5.1 One-soliton solution 



^if-oo \m\ 2 dx' 



2N 



£ — > OO. 



m 



(x)=a^2C je - 2 ^ x M 22 (x,( 3 
3=1 



In this section we derive explicit formulas for the one-soliton solutions. Assume N = 1 
so that there are two zeros of a{Q: Q\ in the first quadrant and £2 = — £i in the third 
quadrant. Using that C\ = C 2 we find that the algebraic system (5.2) reduces to the 
following two equations: 



M 22 (x,Ci) = l + C ie 2 ^ 



1 



1 



Mi 2 (x,Ci) 



,0"Cl-Cl crCl+Cl 
M 12 (x,( 1 ) = C 1 e 2 ^ + -t^-t-) M 22 (x,Ci). 



Solving for M 22 (x, (1) we find 



M 22 (x,Ci) 



A„d\C±\ 2 e 2 ^l-<-> 
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We parametrize the zero Ci an d the normalization constant C\ in terms of four 
parameters A > 0, 7 G (0, it), xq G R, and Xq G K, according to 



(f = A 2 (— CJCOS7 + ishi7), C\ = iA(sin7)e 
Then, writing = A 2 (x — xq) sin 7, we find 



^ifZo \ m ?dx' e 2iSo e 2A 2 a:o sin 7 ^ 



1 + e- 4e e^ CT ^' 



Furthermore, by (5.3), 



m(x) = — 2io"A(sin j)e 4<Tl i*-oo l m l dx 



-2i£+26» 



e 4e + e -<7i 7 • 



where 

Moreover, using 
we deduce that 



= — ctA x cos 7 + So- 

4A 2 sin 2 7 e 4e 

(e 49 + e ai ~y)(e Ae + e~ ai ^) 



2ai[ x \m\ 2 dx' 2<ri[°° \m\ 2 dx> + e ° %1 

O J — CO 1 1 — P J — OO 1 1 

e 40 + eCT i 7 • 



Substituting this into ( |5.3[ ), we find 
(5.4) 

Integration yields 



= 4oA sin 7e 



-2iE+20 e t e 



(TJ7 



U 



( e 40 + eCTi7 )2- 

2isin 7 e - CTi7 e- 2iE 



A e 2e + e CT ^e- 2e ' 
The time evolution is determined by the requirement C\(t) = Ci(0)e 24?? i*: 

S ot = Re [rjl) = -a/3 - aaA 2 cos 7^1 + , 



xot 



Im r/ 2 
A 2 sin 7 



The final expression for the one-soliton solutions of equation (3.1) parametrized by 
the four parameters 



(5.5) 



is 



(5.6) 



7 G (0,tt), 



A > 0, 



Sn G 



x G R, 



2isin7 e -<™7 e -2iE(a;,t) 

A e 29(x,t) _j_ e ai-y e -20(x,t) ' 
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where 



and 



E(./\ / ) — — u SI — rrA' cos -,. ( x + a ( 1 + ) t ) + En 



9(x, t) = A 2 sin 7 ( x — xq + a I 1 



4A 4 



Remark 5.1 The expression for u x as given in equation (5.4) coincides with formula 



(33c) in [TB] for the one-soliton solution of the derivative NLS equation. This is 



because the x-part of the Lax pair (3.2) coincides with the isospectral problem for 



the derivative NLS equation written in terms of u x [16] . 
5.2 Other kinds of soliton solutions 

We finally make some comments regarding the existence of other kinds of soliton 



solutions for equation (1.4). One well-known property of the CH equation is that it 



admits peaked solitons (these are continuous solutions with a peak at their crest) [3]. 



Since (1.4) is related to the NLS equation by a procedure analogous to that which 



gives the Camassa-Holm equation from KdV, it is natural to ask whether (1.4) also 



exhibits weak solutions of this kind. However, a first analysis indicates that there are 
no natural candidates for peakon solutions of equation (1.4) cf. [18 . 



On the other hand, equation (1.4) does admit rational solitons, i.e. smooth soli- 



tons with algebraic decay at infinity. These solutions arise from the solitons with 



exponential decay of equation (|5.6[) in the limit 7 f ir. If u(x,t) is given by (5.6), then 



where 



u r (x, t) 



u(x,t) = U r (x,t) + 0(7T - 7), 7 I 7T, 

2ie~ 2 K A2x+Q (^- ,jA ) 2 ' +<TSo ) A 



4A 4 (x - a? ) - a ((3 2 - 4A 4 ) t + iaA 2 ' 

We claim that u r (x,t) is a rational soliton. Indeed, for the ranges of the parameters 
A, So, xo specified in ( |5.5[ ), u r (x, t) is a smooth function of (x, t) with algebraic decay 
as Ixl 



00, and it can be verified directly that it satisfies equation (3.1). 



A Derivation of a Lax pair 

In this appendix we briefly indicate how the bi-Hamiltonian structure was used to 



derive the Lax pair (3.2). 



1. Consider the recursion operator L = #2#i 1 with adjoint = 1 1 62- The 
bi-Hamiltonian theory implies that the equation 



(A.l) 



JP2 



X 



<P1 

^2 



is the x-part of a Lax pair with spectral parameter A for (1.4) in terms of the 
'squared eigenfunctions' (pi(x,t) and (f2(x,t) cf. [7]. 
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2. Assume the eigenfunctions ipi(x,t) and ^p2(x,t) satisfy an x-part of the form 

f 1pl x = \hli>l + Ui2^2, 
\ ifax = W21"01 + 

for some coefficients hi,h2,ui2,U2i- The 'squared eigenfunctions' (pi = ^1^2 
and (f2 = ^2^1 then satisfy 



flxx + A(/ii - /i 2 ) 



U21x 
U21 



Via 



(A.2) + ( -2ni 2 U2i - (/n - h 2 )\— ) if! + 2w| l¥ > 2 = 0, 



<f2xx + -A(/ii - h 2 ) 



U\2x 



V2a 



U12 

+ ( -2U12U21 + (hi 



U21 



/i 2 )A— ^ ) <p2 + 2ul 2 fi = 0. 

^12 



The coefficients hi, I12, u±2, U21 can now be chosen so that (A.2) coincides with 



(A.l ). After redefining A we find that the x-part takes the form 

I ipix = i{ot\ 2 - f3)ipi + \q1p2, 



(A.3) 



ip2x = -i{a\ 2 - P)ip2 + ArV'i 



for some coefficients a and (5 independent of A. 
3. Assume a £-part of the form 



(A.4) 



■ipit = Aipi + Bi/j 2 , 
ip2t = Ct/ji - Aip 2 , 



and take the coefficients A, B, C to be Laurent series in A. In view of the 



analogous situation for the Camassa-Holm equation (1.3), one might expect 



A, B, and C to include some negative powers of A. After equating terms at each 



order of A, we find that equation (1.4) is the compatibility condition between 



(A.3) and (A.4) provided that the following equations are satisfied: 

ia 1 



A 



and 



a 



X 2 + -£ + Y uv+ ^?T9> B=-q\+— u-, C=-rA+ — v-, 



v 



(77 



P 



1 

2^ 



Av 2 A 2 ' 



q = u + ivu x , 



2v A 



v — IVV X , 



2v A 



aq, v = au. 



4. Transform the Lax pair (A.3)-(A.4) by letting 



(A.5) U 



i' I 

e ix/u v 



. e -ix/v u 




i\fa\, 



V 



e tpi 

ix 

ip2 
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Then equation (1.4) is the condition of compatibility for the two vector eigen- 
value equations 



(A.6) 



V x + i( 2 a 3 V = (U X V, 

Vt + l iU 2 -l + ^) ^ = (%U X + ^a 3 (-U 2 + $U)) V. 



These equations imply the Lax pair (3.2) after suitable redefinitions of a and 
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